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ABSTRACT 

A recent spectrum of the optical afterglow of GRB 970508 suggests that garnma- 
ray bursts (GRBs) are cosmological in origin and it is of crucial importance to derive 
an accurate distance to each burst. If GRBs occur near their host galaxies (<< 40 
kpc) then Lyman limit absorption [IV(HI) > 1.6 x 10 17 cm -2 ] should be observable in 
roughly half the GRB afterglow spectra. Here we outline the methodology to obtain 
a redshift from the GRB afterglow spectrum using the recently installed Space Tele¬ 
scope Imaging Spectrograph (STIS) instrument onboard the Hubble Space Telescope. 
A low-resolution spectrum with the Multi-Anode MicroChannel Array (MAMA) de¬ 
tector gives complete spectral coverage over the wavelength range 1570-3180 A(Near 
UV; NUV) and 1150-1740 A(Far UV; FUV). Assuming a Target of Opportunity obser¬ 
vation is conducted soon (<^ 3 weeks) after a bright burst, a relatively small integration 
time (~ 3 orbits) would be sufficient to detect the Lyman limit over a wide redshift 
range (0.3 <) z <( 2.2). Detection (or non-detection) of the Lyman limit, in concert 
with ground-based observations of nearby galaxies and Mg II and C IV absorption 
lines, should provide meaningful constraints on the relationship of GRBs to galaxies. 
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1 INTRODUCTION 

An optical source associated with gamma-ray burst (GRB) 
970508 has recently been detected (Bond 1997). Absorption 
features seen in a spectrum from the Keck 10m telescope 
indicate that it is at or beyond z = 0.835 and the lack of a 
prominent Lyman-a absorption suggests z < 2.1 (Metzger 
et al. 1997). Presuming that the source is indeed the optical 
afterglow of the burst, then GRBs have only just been con¬ 
firmed to be cosmological. As no emission lines have been 
detected from the transient and it is unclear whether the 
burst is at the redshift of the absorption system at 2 = 0.835 
(see §2), an alternative way to get a limit on the redshift is 
by looking for a Lyman -a forest or Lyman limit absorption. 
Lyman limit absorption arises from neutral hydrogen (HI) 
which is optically thick to Lyman-continuum radiation for 
wavelengths A < 912 A in the rest frame of the absorbing 
system. Since both redshifted Lyman-a absorption (1216 A) 
and the Lyman limit (A < 912 A) remain in the UV pass- 
band for 2 < 2.2, currently only the Space Telescope Imaging 
Spectrograph (STIS) onboard the Hubble Space Telescope 
(HST) can detect a Lyman limit. As expected from theoret¬ 
ical models, both optical transients detected thus far have 
started fading soon after discovery, and thus a STIS spec¬ 


trum must be obtained as soon as possible after the burst 
so as to maximise signal-to-noise (S/N). 

The distance inferred from a Lyman limit in the con¬ 
tinuum would provide knowledge of the luminosities and en¬ 
ergies involved in the explosion—two vital parameters for 
constraining GRB models. In section 2 we discuss the pos¬ 
sibilities of the GRB afterglow undergoing Lyman limit ab¬ 
sorption. Then, in section 3 we discuss the instrumentation 
and calculate the integration time required to infer the red¬ 
shift of a GRB. 


2 LYMAN LIMIT ABSORPTION 
EXPECTATIONS 

Lyman limit absorption systems, which are generally HI 
clouds opaque (r > 1) to the Lyman-continuum, are be¬ 
lieved to concentrate in the disk and halo of most galaxies. 
Steidel (1993) has found that the density of Lyman limit sys¬ 
tems is high enough out to impact parameter of ~ 40 kpc in 
the galactic disk that continuum radiation passing through 
the plane of the disk will always be subjected to Lyman 
limit absorption. Thus, if GRBs occur at offsets <<40 kpc 
from the centres of their host galaxy, it is expected that 
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roughly 50 percent (half in front of the disk, half behind the 
disk) of GRB afterglow will have a Lyman limit break in 
the spectrum; this Lyman limit will correspond to the pre¬ 
cise redshift of the burst since the limit system will be local 
to the GRB. 

Intervening galaxies, not associated with the GRB but 
in the line-of-sight of the afterglow, may also absorb the 
continuum; thus a redshift inferred from a Lyman limit will 
not necessarily be the redshift of the GRB. What effect will 
this have on the determination of GRB redshifts? Storrie- 
Lombardi et al. (1994) survey QSO absorption spectra and 
find that for any random line-of-sight, the density of inter¬ 
vening Lyman limit systems is N(z) ~ 0.38(1 + z) 1 ’ 04 for 
redshifts z < 3.0. Thus at redshifts z > 1.5 it is expected 
that most GRB afterglows will be subjected to at least one 
Lyman limit in their continuum that does not necessarily 
correspond to the intrinsic redshift of the GRB. 

If GRBs are ejected to distances comparable to the scale 
length of the Lyman limit absorption systems in the disk, 
then the probability that the host galaxy will absorb the 
spectrum shortward of the Lyman limit is reduced. The top 
portion of figure (jl:) shows the expected probability of the 
existence of a Lyman limit (r > 1) in the spectrum of a GRB 
afterglow as a function of redshift and the offset scale length 
of GRBs from their host galaxy. The relationship between 
the frequency of absorption from the host galaxy and the 
offset scale is computed by assuming that absorption only 
occurs if the GRB is seen through the 40 kpc absorbing disk 
and that the disk has random viewing inclination. As seen, 
the frequency of Lyman limit absorption in the spectrum of 
GRBs at low redshifts (z < 1) may be used to determine 
the intrinsic offset of GRBs from their host galaxies since 
most Lyman limit absorption at low redshifts is expected to 
come the host galaxy. The bottom half of figure (jlj) shows, 
as a function of offset scale and redshift, the probability that 
the inferred redshift is within 20 percent of the redshift of 
the GRB afterglow. If GRBs occur with about 60 kpc of 
their host galaxy (solid line), then more than 60 percent of 
the redshifts inferred from the Lyman limit in the spectrum 
will be a moderately accurate (< 20 percent) measure of the 
redshift of the GRB. 

Absorption from Mg II, with a galactic impact param¬ 
eter of ~ 50 kpc (Bergeron et al. 1994), and C IV is ex¬ 
pected, but not required, to accompany a Lyman limit sys¬ 
tem (Storrie-Lombardi et al. 1994). Thus, figure ( lj) could 
also be seen as a prediction of the frequency of absorption 
lines in GRB afterglows. Indeed for GRB 970508, both Mg 
II and Fe II absorption was detected using Keck (Metzger et 
al. 1997); Arav & Hogg (1997) have found that the absence 
of detectable C IV in the spectrum of the afterglow limits 
the redshift of the GRB to z < 1.8 (95 percent confidence). 

3 OBSERVING DETAILS 

The STIS CCD instrument onboard HST with a G230LB 
low-resolution grating gives complete spectral coverage over 
the wavelength range 1685-3065 A. This would be suffi¬ 
cient to detect the Lyman limit over the redshift range 
(0.85 < z < 2.4). The preferred instrument for obtain¬ 
ing redshifts from UV spectra, however, are the Near UV 
(NUV) and Far UV (FUV) Multi-Anode MicroChannel Ar- 
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Figure 1. Top plot: The expected probability of the existence a 
Lyman limit in the spectrum of a GRB afterglow as a function of 
redshift and the offset scale length of GRBs from their host galaxy. 
The contribution to the total probability of a Lyman limit system 
(LLS) from the host galaxy (0.5 solid; 0.3 dashed; 0.1 dot dashed) 
is computed by assuming the GRBs occur randomly around the 
galaxy and the impact parameter of the LLSs in the disk is ~ 40 
kpc (Steidel 1993). The contribution from intervening LLS (not 
associated with the GRB itself) is calculated assuming the num¬ 
ber density of LLS systems evolve as N(z) = 0.38(1 + z ) 1,04 
(Storrie-Lombardi et al. 1994). Bottom plot: The probability, if a 
Lyman limit is present in the spectrum of the afterglow, that the 
inferred redshift is accurate to within 20 percent of the true red¬ 
shift of the GRB afterglow itself. Not surprisingly at low redshift, 
a distance inferred from the Lyman limit is probably the true 
distance to the GRB since the expected number of LLSs from 
intervening galaxies is small. See Appendix A for details. 

ray (MAMA) detectors which is more sensitive and has much 
wider spectral range (FUV: 1150-1736 A; NUV: 1570-3180 
A) corresponding to a broad redshift range (0.26 < z < .9 
and 0.72 <) z <( 2.5, respectively). 

3.1 Signal—to—Noise (S/N) estimation 

Temporal fits to the broadband spectra of the optical tran¬ 
sients GRB 970228 (Van Paradijs et al. 1997) and GRB 
970508 (Bond 1997) indicate that the GRB afterglow may 
follow a power-law both in flux and spectral shape in optical 
wavelengths. It should be noted that the later transient had 
an observed rising phase of optical emission (e.g. Galama et 
al. 1997). Although individual cases may vary, in calculating 
the expected signal as a function of wavelength and time, we 
adopt the functional form of the flux as: 

F v (t) = F to u~ s (t - to)' 3 (1) 

where Ft 0 properly normalises the spectrum at to, the time 
at which the decay begins. Wijers, Rees, & Meszaros (1997) 
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Figure 2. The expected integrating time required to achieve a 
S/N=3 in a 100 A bin as a function of U magnitude of the after¬ 
glow and redshift of the Lyman limit source. The MAMA detec¬ 
tors on STIS (NUV, solid lines; FUV dashed) achieves a signif¬ 
icantly better signal-to—noise than the CCD (dotted lines). For 
both MAMA NUV and the CCD, the redshifts correspond to (left 
to right): z = 0.9, z = 1.2, z = 1.5, and z = 2.1. The FUV curves 
correspond to (left to right): z = 0.3, 0 = 0.7, and z = 0.5. Also 
shown on the right vertical axis is the number of HST orbits re¬ 
quired assuming the object is observable for 90 minutes of each 
orbit. The signal—to-noise is calculated for the indicated grating 
with a 52x 0.2 arcsec slit and assuming average Zodiacal light, 
bright earth level. 


find that the afterglow adequately fits the early light curve 
with 5 = 0.8 and /3 = —1.2 for GRB 970228. Our preliminary 
fits to the data from GRB 970508 indicate that 5 ~ 0.8 
and /3 ~ —1.0. If the afterglow is observed as a Target of 
Opportunity t 0 b s < 3 weeks after the burst (e.g. U(t 0 b s ) ~ 
22.7 for GRB 970508), a redshift could be obtained in less 
than 10 HST orbits (see fig. [2]) using the either FUV or 
NUV MAMA detectors. 

Unless HST is able to observe the afterglow of a GRB 
while it is still bright (U <j 20), detection of damped Ly¬ 
man a absorption at low redshift (z < 1.5) will require a 
very large integration time on STIS. However, a significant 
detection of a Lyman limit requires far less S/N per unit 
wavelength and thus improves the chance of determining a 
redshift of GRBs from fewer orbits. Although the Lyman 
break occurs at shorter wavelengths (A < 912(1 + z) A) 
than the Lyman a forest (A ~ 1216(1 + z) A) where the de¬ 
tectors are less sensitive, the distinct cut-off of this spectral 
feature is unambiguous and does not require good spectral 
resolution, making this the most effective and clearcut way 
to limit the redshift of faint sources. 

Figure (|i|) shows the expected integration time required 
to achieve a S/N=3 in a 100 A bin as a function of U magni¬ 
tude of the afterglow and redshift of the Lyman limit source 


Figure 3. Simulated afterglow spectrum with ~ 1.5 HST orbits 
(5400 sec). The spectrum is computed for the MAMA 52x 0.2 
arcsec G230L slit on STIS; the source is taken to be a magnitude 
U = 21.0 and the spectral shape is <5 = 1 (eq. 1). We have re¬ 
moved from the spectrum the narrow Oil line (2470 A) expected 
from geocoronal emission. The Lyman limit is clearly detected at 
A ~ 2100 A corresponding to a redshift z — 1.3. The histogram 
overlay, with estimated error bars, is the renormalised counts per 
100 A with arbitrary offset. 


for both the CCD and MAMA detector on STIS^J As seen 
in the figure, if the afterglow is observed while it is still 
reasonably bright (U < 21), detection of a redshift z > 0.3 
would require less than 1 HST orbit. Figure (||) shows a sim¬ 
ulated spectrum of an afterglow obtained with MAMA with 
~ 1.5 HST orbits (5400 sec) where the source is a magni¬ 
tude U = 21.0 and the spectral shape is 5 = 1 (eq. 1). We 
have removed the contribution of the source to the observed 
spectrum for wavelengths A < 2098 A corresponding to a 
Lyman limit at a maximum redshift of z = 1.3. The clear 
break in the spectrum implies that a fairly accurate redshift 
is obtainable. 


4 DISCUSSION 

In general, given the efficiency of STIS and the spectral 
shape of the afterglow, the S/N increases with higher red¬ 
shift (see fig. Q) for the NUV MAMA detector. If the burst 
occurs at a redshift z )> 1.5 then, with only a few orbits, 
one could determine the position of the Lyman break well 
enough to determine its redshift to an accuracy of better 
than 5 percent; it may even be possible to detect a Lyman- 
a forest at high z with only a few orbits. As it is easier to 


* Calculated usin; 
culator located at 


the STIS Sncctrosconic Exposure Time Cal- 
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detect a redshift for bursts that originate from a higher z 
with the NUV MAMA detector, the absence of a Lyman 
limit in the spectrum would, in general, always provide an 
upper limit to the redshift. From figure (Q), it is clear that 
use of the MAMA detectors onboard the HST, given their 
very low internal countrates, are preferable to the CCD. 

The disadvantage of using HST to infer a redshift is the 
difficulty of altering its scheduled observations even on the 
timescale of weeks. Certainly the advantage is that there are 
very few available ground-based telescopes that can resolve 
Mg II and C IV absorption lines from faint objects; even 
then, such telescopes may not have the proper viewing con¬ 
ditions to the source and detected lines may place different 
limits in the redshift. In addition, there may be important 
emission (e.g. lines) in the UV spectrum of GRBs that would 
not be observable from the ground. 

As the ensemble of optical counterparts begins to grow 
in size, it will be possible to infer both their redshift distri¬ 
bution and the distribution of GRBs with respect to observ¬ 
able galaxies by noting the frequency of absorption lines and 
detected Lyman limit (see fig. [I]). If no Lyman limit is de¬ 
tected in the spectra of the optical transients (especially at 
low-redshifts), the conclusion would be that either GRBs do 
not originate near galaxies (< 300) kpc or they are Galactic 
in origin. 


5 CONCLUSIONS 

A precise redshift of a gamma-ray burst would greatly fur¬ 
ther the field by providing an accurate understanding of the 
energies involved in the burst. Assuming the extrapolation 
of the spectral index (5) from the optical to UV passband is 
correct, we find that a MAMA detector observation of the 
bright afterglow of a GRB by STIS over ~ 1 orbit could re¬ 
veal a Lyman limit in the spectrum and hence, if the GRB 
occurs behind a region of neutral hydrogen of its host galaxy, 
provide a direct redshift to a gamma-ray burst. Given that 
Lyman limit absorption may come from galaxies along the 
line-of-sight, we predict the expected frequency of Lyman 
limits in the spectrum as a function of GRB redshift and 
the distance from their host galaxy. 
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APPENDIX A: LLS PROBABILITY 

A Lyman limit system (LLS) is defined as the region of space 
for which the optical depth (r) to the Lyman-continuum is 
greater than one. Most Lyman limit systems are believed to 
be clouds of HI with a column density 7V(HI) > 1.6 x 10 1 ' 
cm 2 . As depicted in the top half of figure 1, the probability 
as a function of redshift that a spectrum of a GRB will 
contain at least one such LLS is 

P(> 1 LLS) = 1 — Poisson [0 LLS | m(z m i n , z)] 

x P(No LLS from host Galaxy) 

where Poisson [0 LLS|m(z m in, a)] is the Poisson probability 
of no intervening LLS between the GRB source (at redshift 
z) and the observer given the expected number of LLSs: 

m(z m in,z)= f N(z')dz'. (2) 

“ 2 min 

The minimum redshift in which a LLS could be detected 
is 2min and N(z) is the number density of LLSs per unit 
redshift. 

The probability that the redshift inferred from an ob¬ 
served Lyman limit in a GRB afterglow spectrum is at least 
x times the redshift of the GRB itself is given as 

P(^lls > xxz \>1 LLS) = where, (3) 

P(A) 

P(A | B) = P(> 1 LLS| LLS redshift > x x z) = 1, (4) 

P(A) = P(> 1 LLS), (5) 

P(B) = P(z L ls > x x z) (6) 

= 1 — Poission [0 LLS | m(x x z, z)] 

x P(No LLS from host Galaxy). 

This probability is depicted in fig. 1 for x = 0.8. 
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